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Acoustic Energy Storage in Single Bubble Sonoluminescence 
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Single bubble sonoluminescence is understood in terms of a shock focusing towards the bubble 
center. We present a mechanism for significantly enhancing the effect of shock focusing, arising from 
the storage of energy in the acoustic modes of the gas. The modes with strongest coupling are not 
spherically symmetric. The storage of acoustic energy gives a framework for understanding how light 
intensities depend so strongly on ambient gases and liquids and suggests that the light intensities of 
successive flashes are highly correlated. [S0031-9007(96)01412-3] 



PACS numbers: 78.60.Mq, 42.65.Re, 43.25. + y, 47.40.Nm 

Sonoluminescence (SL), the conversion of acoustic en- 
ergy into light, occurs by coupling gaseous bubbles to an 
externally forced liquid [1]. Two experimental configura- 
tions for SL exist: multibubble sonoluminescence (MBSL) 
[1,2], which occurs in transient cavitation clouds, and sin- 
gle bubble sonoluminescence (SBSL), which occurs when 
a single bubble is trapped at the node of an applied acoustic 
field [3,4]. Recent experiments [3-10] uncovered many 
remarkable properties of SBSL, including picosecond light 
pulses [4,11] and sensitive dependences on almost all ex- 
perimental parameters. 

SL requires both energy transfer from the liquid to the 
gas, and focusing of this energy. The shock theory stipu- 
lates that a shock focuses the energy input during a single 
collapse. Large temperatures occur because a focusing 
shock is a singular solution to the Euler equations [12] in 
which the maximum temperature diverges at the focusing 
point. Although this shock theory gives consistent expla- 
nations for many aspects of SL (e.g., high temperatures and 
picosecond light pulses), there is mounting experimental 
evidence that tensions in the theory exist. To wit: While 
the shock theory is believed to work for both MBSL [13] 
and SBSL [14], the entire MBSL bubble cloud emits less 
than 1 % of the light of a single bubble in SBSL. For SBSL 
in alcohols, Weninger et al. [9] reported the existence of 
abrupt jumps in the light intensity, in which the light output 
changes by a factor of 400 with a 1 °C change in the liquid 
temperature. More recently, Weninger et al. [10] reported 
the existence of angular correlations in the light output, in 
which the radiation field had a significant dipole moment. 
We believe that it is difficult to resolve these discrepancies 
within the single shock theory, and that additional physics 
is needed. 

This paper presents a mechanism for enhancing the ef- 
fect of shock focusing, giving natural resolutions to the 
aforementioned tensions in the shock theory. The idea is 
that the energy focused in SBSL is not input in a single 
bubble collapse, but instead accumulates within the acous- 
tic modes of the gas over many bubble cycles. Because 
the energy stored in the modes is far greater than the en- 
ergy input in a single collapse, the focusing power of the 



shocks is significantly enhanced. (Note that the emitted 
light carries away very little energy, and merely acts as an 
indicator of the energy density at collapse [14].) Within 
this picture, the maximum temperature and light intensity 
of SBSL are largely set by the total acoustic energy stored 
within the bubble, instead of the details of the focusing sin- 
gularity. An essential condition for accumulative energy 
storage is the stability of the bubble [4] from surface, dif- 
fusive, and chemical instabilities [6,15,16], which exists in 
SBSL but not in MBSL. 

Our theory is motivated by the fact that the basic 
assumption of the shock theory — that the shock launches 
into a gas with a spatially uniform density — breaks down 
in the absence of dissipation. Without dissipation, energy 
transferred from the liquid to the gas in a single collapse 
never leaves the bubble, resulting in a nonuniform gas 
density over long times. This inconsistency is cured 
only with enough dissipation so that the gas dynamics is 
overdamped. When the system is underdamped energy 
builds up in the acoustic modes of the bubble. 

The key quantity distinguishing the underdamped and 
overdamped regimes is the Floquet multiplier A, char- 
acterizing the net input of energy to the acoustic modes 
in a single cycle. To be self-consistent, the shock the- 
ory requires A < 1; acoustic energy storage corresponds 
to A > 1. The total multiplier is the product A = AD 
of multipliers for energy input A and dissipation D < 1. 
Energy input occurs mainly at the bubble collapse through 
resonances and shocks. Dissipation occurs throughout the 
cycle and is the product of several factors, mainly viscous 
dissipation in the gas and acoustic energy transmission 
from the bubble to the fluid. 

In the following, we first compute the acoustic modes 
of the gas, and then address A for the different modes. It 
emerges that the two most easily excited modes are not 
spherically symmetric. Then, we discuss dissipation. We 
argue that varying the dissipation naturally leads to tran- 
sitions between the underdamped A > 1 and overdamped 
A < 1 regimes, manifest by an abrupt jump in the light 
intensity, as seen by Weninger [9]. Finally, we present 
further comparisons and predictions for experiments. 
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We now turn to a calculation of the normal modes of the 
gas. To keep it simple, assume nondissipative, isothermal 
gas dynamics in the bubble with a velocity v(r, t) and gas 
density p(r, t), 



d t p = —V • (pi;), d t v + v • Vv = —c 
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gas 
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Here, c gas is the (constant) sound speed which is deter- 
mined by the gas temperature & g . The boundary condi- 
tion for the velocity is set by the Rayleigh-Plesset equation 
[6,17] 
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which describes the dynamics of the bubble wall radius 
Rit), with viR,t) ■ r = R(t). Typical parameters for 
an argon bubble in water [8] are surface tension a = 
0.073 kg/s 2 , water viscosity v = 10~ 6 m 2 /s, density 
pi = 1000 kg/m 3 , and speed of sound c/ = 1481 m/s. 
The driving frequency of the acoustic field Pit) = 
P a costot is cj/2tt = 26.5 kHz [8] and the external 
pressure Pq = 1 atm. The pressure inside the bubble 
varies as piR) « (7? 3 — h 3 )~ K , with h the hard core van 
der Waals radius. Here k is the effective polytropic 
constant, determined by thermal effects to be k ~ 1 [18]. 

The density and velocity consistent with the Rayleigh- 
Plesset equation are 
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where M is the Mach number M = R/c gas . The normal 
modes follow from linearizing about this state by writing 
p = p s + p(r/R,t) and v = v s + i>ir/R,t). Then 
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where x = r/R, and V x is the gradient with respect to x. 
Introducing the velocity potential via vix,t) = VO(r, t) 
yields a wave equation for t), 



d 2 = c 



2 yl 
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The general solution to this equation is a sum of modes 
®(r, t) = ai,n(t)ji(ki t „r/R)Yi, where j t are spherical 
Bessel functions, 7/ are spherical harmonics depending 
only on the angle variables, and the „ are determined by 
the boundary condition for the radial velocity i>iR,t)f = 
0, i.e., the rath zeros = 0. The coefficients a/_„ 

satisfy the equation for a linear oscillator with time 
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varying frequency (often called Hill equation) 

+ ^l,n a hn =0, (7) 

with ft />B = c gdLS ki t „/R. 

When there is a large separation of scales between 
the oscillator frequency „ and R/R, no energy is fed 
into a/ „; the oscillator energy divided by its frequency, 
and therefore the envelope / of a 2 n /R is an adiabatic 
invariant [19]. When the separation of scales is violated 
(R ~ Cgas^/n), energy is transferred into a; „ and / is 
no longer invariant, but increases by a factor A. The 
most efficient energy transfer (largest A) happens for the 
modes with the smallest values of ki >n . The first three are 
Jki.i « 2.0816, k 2 ,i « 3.3421, and %i « 4.4934. Thus, 
the two most easily excitable modes are not spherically 
symmetric [20]. 

To demonstrate the energy buildup in acoustic modes, 
Fig. 1 shows the most unstable mode a\ t \{f), together with 
Rit) and a\\/R, clearly demonstrating the parametric 
instability of Eq. (7). The fast time scale in aij(f) 
corresponds to the acoustic travel time. 

This energy buildup eventually saturates due to dis- 
sipation and nonlinear effects. We do not understand 
the detailed structure of the saturated state. However, in 
general, many different modes a/ „ will be excited and in- 
teract. The shock produced by the bubble collapse will fo- 
cus the total energy contained in this saturated state. This 
shock is much stronger than when only the energy input 
during a single collapse is focused; thus, the light intensity 
is significantly enhanced. Since shock focusing in spheri- 
cal geometry is genetically singular, we still expect that 
there are qualitative similarities between shock focusing in 
this highly nonuniform gas and shock focusing in a homo- 
geneous gas [12,14]. In particular, the width of the light 
pulse within this picture should be comparable to the width 
in the shock picture. 
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FIG. 1. (a) Solution Rit) to the Rayleigh-Plesset equation. 
The forcing pressure is P a = 1.4 atm and the ambient radius 
Ro = 7 pm. (b) Buildup of energy in the normal mode ai,i(0- 
(c) aj y i/R whose envelope is the adiabatic invariant /. 
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We now proceed to analyze parameter dependences 
of this acoustic resonator mechanism. The above model 
neglects dissipation. However, in the real system two 
major types of dissipation are present, namely, viscous 
dissipation and acoustic wave transmission from the gas 
into the liquid. First we consider viscous dissipation: 
a viscosity /jl causes the adiabatic invariant to decrease 
between successive bubble collapses by a factor D v . For 
a mode of wave number k, I puj p(k/R) 2 I , so that 



D v ~ exp 



/„ 



/* k 2 



(8) 



where T is the bubble oscillation period. Since the bubble 
spends most of its time near the maximum radius R max , 
we estimate D v ~ exp[— T£ 2 /ii//?^ ax p(/? max )]. A crucial 
point is that p(R max ) is not determined by Eq. (3) for 
p s but instead by the vapor pressure p v of the liquid. 
This parameter depends on the liquid temperature 0/ via 
p v °c e -Ejk B @, ( w i t h £ a an activation energy). 

Another dependence of the Floquet multiplier A on 0/ 
comes through A. The energy input A becomes larger for 
smaller c gas , since decreasing c gas increases the time that 
the liquid and gas frequencies are comparable. Numerical 
simulations of Eq. (7) show that A decays exponentially 
with c gas . Since the gas temperature before the bubble 
collapse is equal to the liquid temperature [6,18], we 
have c gas K V®s = Combining these results, we 

can now plot the 0; dependence of the total Floquet 
multiplier A(0 ; ) = A(0/)D„(0/), which is shown in 
Fig. 2(a) [21]. When 0/ is decreased, A increases up to 
some temperature where the light intensity is maximal. 
At a critical temperature 0;* [indicated by the arrow 
in Fig. 2(a)] dissipation beats the input and a transition 
towards the overdamped regime occurs, manifested by an 
abrupt change in the acoustic energy stored by the bubble, 
and thus in the light intensity. 

As mentioned above, this kind of behavior is found by 
Weninger et al. [9] in experiments on alcohols [Fig. 2(b)]. 
Decreasing 0/, the light intensity increases to a maximum, 
followed by an abrupt transition. In this transition, the 
light intensity decreases by a factor of 400 with a 1 °C 
change in 0/. The above argument suggests that the 
crossover temperature 0;* should be largest for the fluid 
with the smallest vapor pressure, all other parameters 
being constant. This is confirmed by the experiments [9]: 
of the five alcohols tried, transitions occur in 1-pentanol 
and 1-butanol, with the lowest vapor pressures at room 
temperature (0.003 and 0.007 atm, respectively) [22]. 

Slightly above the temperature ©;* of the abrupt tran- 
sition, the Floquet multiplier is A ~ 1 + a(&i — ©/). 
The number of cycles N it takes for the system to satu- 
rate satisfies A N £ sat , with £ sat the saturation energy. 
This gives that N diverges like N (0; - 0;*) _1 . The 
theory also suggests that below the sharp transition the 
SBSL intensity should be comparable to that in MBSL, 
since single shocks are responsible for the light intensity. 
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FIG. 2. (a) A sketch of the qualitative behavior of the Floquet 
multiplier A(® ; ). Scales are adopted to match (b): Depen- 
dence of SL intensity on @ ; for a xenon gas bubble at a partial 
pressure of 150 mm Hg for two alcohols, as measured by [9]. 
The transition turning off the collective acoustic state is marked 
on the graph. 1-propanol shows no transition in this tempera- 
ture regime because its vapor pressure is larger than that of 
1-pentanol. Similar transitions have also been observed in bu- 
tanol at slightly different parameters. 



We also speculate that, due to the smaller central tempera- 
tures of bubbles below the transition, SBSL spectra may 
show more structure in this regime, perhaps resembling 
those of MBSL. 

Whether or not the abrupt transition actually occurs 
in a given system depends on (a) the value of A at the 
liquid freezing-boiling point (i.e., whether the transition 
is in the liquid regime of the fluid), and (b) whether or 
not long time bubble stability is affected by changing 0/ 
(e.g., by changes in the gas solubility in the liquid which 
is a central parameter for stability [15]). The fact that 
this transition happens in the alcohols demonstrates that 
energy input to the gas by the Rayleigh-Plesset dynamics 
is sufficiently high for A > 1 to occur. Because bubbles 
in water are much brighter than bubbles in alcohol, we 
conjecture that this common situation is also underdamped. 
One may experimentally test the conjecture by carefully 
examining the temperature dependence of the onset of 
SBSL in air bubbles in water: At a critical driving pressure 
P* (depending on 0;, c x , etc.), the light turns on. By 
keeping the driving pressure fixed slightly above P* and 
then varying the liquid temperature 0;, the light should 
turn off abruptly, in a fashion similar to Fig. 2. 

There are other contributions to A besides those men- 
tioned above. An important loss mechanism is acoustic ra- 
diation through the wall, which leads to a reduction of the 
stored energy by a factor D a ~ I — (p s c gas )/(p/C;) for 
every wave bounce in the bubble. When the bubble is fully 
expanded, D a ~ 1 since p s <C pi and no losses occur. 
However, significant transmission occurs when the bubble 
is maximally compressed, where the gas density is of order 
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the liquid density. Therefore, the light intensity should in- 
crease with the density of the liquid, and decrease with the 
density of the gas (with all other parameters held constant). 
Both these effects have been observed: experiments in al- 
cohols [9] show a systematic increase in the light intensity 
with increasing density of the liquid. At room temperature 
the heaviest alcohol, pentanol, emits about 100 times more 
light than the lightest alcohol, ethanol. As for the gas in 
the bubble, the spectrum of xenon bubbles has a peak in 
the visible, whereas the lighter helium and argon bubbles 
show no such peak [5], suggesting a larger total intensity 
for the lighter gases. 

Another experiment naturally interpreted is Weninger 
et al.'s observation [10] that the light emission from the 
bubble can be dipolar, not spherically symmetric. The 
dipole moment of the light emission shows significant 
correlations for delay times as long as 100 bubble os- 
cillation cycles, many orders of magnitude longer than 
the pulse width. The present theory predicts that the gas 
density should have a significant dipole component, since 
the most easily excited mode shows this symmetry. This 
nonuniformity in the density causes differences in the in- 
dex of refraction, producing an anisotropic light emission. 
The strength of the anisotropy depends on the size of the 
nonspherical component when the parametric instability 
saturates. 

A basic difference between the present theory and the 
shock theory is that whereas the shock theory provides 
that successive flashes of light are completely independent, 
the present theory predicts strong correlations. Within 
our interpretation, the long lifetime of the dipole state 
measured by Weninger et al. [10] reflects these long time 
correlations. 

To conclude, we have presented a mechanism for en- 
hancement of shock focusing, and argued that this en- 
hancement occurs in SBSL. The mechanism requires long 
time bubble stability, and thus does not work for transient 
bubble collapse in MBSL. By varying dissipation, the en- 
hancement mechanism can be turned off, thus leading to 
abrupt jumps in light intensity [9]. The nonspherical gas 
motion implies a natural explanation for a long lived dipo- 
lar component of the light emission [10]. 

We thank S. Grossmann and L. Levitov for helpful 
discussions, and L. P. Kadanoff for comments on the 
manuscript. M. B. acknowledges an NSF postdoctoral 
fellowship; R. R. R. was partially supported by the NSF 
DMS Grant No. 9311438; M.B. and R.R.R. acknowl- 
edge the Sloan Fund of the School of Science at MIT. 
S. H. and D. L. acknowledge support by the DFG through 
itsSFB185. 



[1] S. Putterman, Sci. Am. 272, 32 (1995); L. A. Crum, Phys. 

Today 47, No. 9, 22 (1994). 
[2] H. Frenzel and H. Schultes, Z. Phys. Chem. 27B, 421 

(1934). 



[3] D.F. Gaitan, L.A. Crum, R. A. Roy, and C.C. Church, 

J. Acoust. Soc. Am. 91, 3166 (1992). 
[4] B.P. Barber and S.J. Putterman, Nature (London) 352, 

318 (1991). 

[5] R. Hiller, S.J. Putterman, and B.P. Barber, Phys. Rev. 
Lett. 69, 1182 (1992); R. Hiller, K. Weninger, S.J. 
Putterman, and B.P. Barber, Science 266, 248 (1994); 
R. Hiller and S.J. Putterman, Phys. Rev. Lett. 75, 3549 
(1995). 

[6] R. Lofstedt, B.P. Barber, and S.J. Putterman, Phys. Fluids 
A 5, 2911 (1993); B.P. Barber et al., Phys. Rev. Lett. 72, 
1380 (1994). 

[7] B.P. Barber, C.C. Wu, R. Lofstedt, P.H. Roberts, and 

S.J. Putterman, Phys. Rev. Lett. 72, 1380 (1994). 
[8] B.P. Barber, K. Weninger, R. Lofstedt, and S.J. Put- 
terman, Phys. Rev. Lett. 74, 5276 (1995); R. Lofstedt, 
K. Weninger, S.J. Putterman, and B.P. Barber, Phys. Rev. 
E 51, 4400 (1995). 
[9] K. Weninger et al, J. Phys. Chem. 99, 14 195 (1995). 

[10] K. Weninger, S. Putterman, and B. Barber, Phys. Rev. E 
(to be published). 

[11] M.J. Moran et al, Nucl. Instrum. Methods Phys. Res., 
Sect. B 96, 651 (1995). 

[12] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perga- 
mon Press, Oxford, 1987). 

[13] P. Jarman, J. Acoust. Soc. Am. 32, 1459 (1960). 

[14] H. P. Greenspan and A. Nadim, Phys. Fluids A 5, 1065 
(1993); C.C. Wu and P.H. Roberts, Phys. Rev. Lett. 70, 
3424 (1993); Proc. R. Soc. London A 445, 323 (1994); J.I. 
Gersten, L. Kondic, and C. Yuan, Phys. Rev. E 52, 4976 
(1995); W. Moss, D. Clarke, J. White, and D. Young, 
Phys. Fluids 6, 2979 (1994); V.Q. Vuong and A.J. Szeri 
(to be published). 

[15] S. Hilgenfeldt, D. Lohse, and M.P. Brenner, "Phase 
Diagrams for Sonoluminescing Bubbles," Phys. Fluids 
(to be published); M. Brenner, D. Lohse, D. Oxtoby, 
and T. Dupont, Phys. Rev. Lett. 76, 1158 (1996); M.P. 
Brenner, D. Lohse, and T. F. Dupont, Phys. Rev. Lett. 75, 
954 (1995). 

[16] D. Lohse et al., "Single Sonoluminescing Bubbles as 

Stable Chemical Reaction Chambers" (to be published); 

M. P. Brenner, S. Hilgenfeldt, and D. Lohse, in Nonlinear 

Physics of Complex Systems — Current Status and Future 

Trends, edited by J. Parisi, S. C. Miiller, and W. Zimmer- 

mann (Springer, Berlin, 1996). 
[17] L. Rayleigh, Philos. Mag. 34, 94 (1917); M. Plesset and 

A. Prosperetti, Annu. Rev. Fluid Mech. 9, 145 (1977); 

W. Lauterborn, J. Acoust. Soc. Am. 59, 283 (1976). 
[18] V. Kamath, A. Prosperetti, and F. Egolfopoulos, J. Acoust. 

Soc. Am. 94, 248 (1993). 
[19] E.J. Hinch, Perturbation Methods (Cambridge University 

Press, Cambridge, 1991). 
[20] This observation, in itself, invalidates the spherically 

symmetric shock theory usually assumed. 
[21] There may be other effects contributing to the ©; 

dependence. 

[22] A nonalcoholic fluid (n-dodecane) has a lower p v and 
shows no transition. No data are available below °C 
because of bubble instability. 



3470 



